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I.  INTRODUCTION 


A.  Study  Objective 

The  purpose  of  this  study  is  to  investigate  the  poten¬ 
tial  existence  of  a  clustering  of  system  poles  and  zeros  in 
the  z-domain.  If  the  potential  exists,  a  mathematical 
means  of  obviating  the  problem  is  to  be  determined. 

The  study  effort  is  composed  of  a  single  task  described 
in  the  Scope  of  Work  contained  in  Technical  Requirement  No. 
T-0208  entitled  "Application  of  Digital  Technology  to  Guid¬ 
ance  and  Control  Theory"  dated  21  May  1979  [1] .  The 
requirements  imposed  by  this  task  are  contained  in  the 
referenced  document  and  re-stated  (for  the  convenience  of 
the  reader)  below. 

Task  I  requirements: 

2.1  Investigate  to  determine  if  a  sensitivity 
problem  due  to  clustering  of  the  system  poles  and 
zeros  exists  for  the  parameter  space  technique  of 
digital  system  analysis/design. 

2.2  If  a  sensitivity  problem  (i.e.,  cluster¬ 
ing  of  the  system  poles  and  zeros)  is  found  in  the 
investigation  of  Paragraph  2.1,  a  technique  shall 
be  developed,  to  mathematically  spread  apart  the 
system  poles  and  zeros  via  intermediate  complex 
transform  techniques  to  eliminate  the  sensitivity. 

B.  Study  Schedule 

The  requirements  of  Task  I  were  completed  during  the 
period  1  August  1979  (contract  initiation)  through  30 
September  1979.  They  are  reported  upon  in  this  report. 

C.  Contents  of  Report 

This  technical  report  covers  the  efforts  performed  in 
completing  the  requirements  of  Task  I  of  the  study.  As 
stated  in  the  referenced  contract  requirements:  "The 
required  documentation  shall  be  a  technical  report." 
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II.  PROBLEM  INVESTIGATION 


As  indicated  in  Paragraph  2.1  of  the  referenced  Tech¬ 
nical  Requirement,  an  investigation  was  made  to  determine 
if  a  sensitivity  problem  due  to  clustering  of  the  system 
poles  and  zeros  exists  for  the  parameter  space  technique  of 
digital  system  analysis/design.  As  will  be  shown  in  the 
sequel,  any  technique  that  utilizes  design  in  the  complex 
z-domain  is  subject  to  this  problem.  If  this  problem  is 
shown  to  exist,  then  the  parameter  plane  technique  (among 
others)  is  subject  to  it  [2] .  Hence  it  is  prudent  to 
describe  what  the  nature  of  the  problem  is  (Section  A, 
below)  and  see  how  likely  it  is  to  occur  in  practice 
(Section  B,  below). 

A.  Nature  of  Problem 

The  nature  of  the  problem  can  be  seen  in  the  complex 
z-domain.  Consider  that  the  dynamics  of  a  system  are 
represented  by  three  roots,  such  as  system  poles:  two  that 
are  complex  conjugates  and  one  that  is  real.  Further 
suppose  that  they  are  all  located  arbitrarily  close  to  the 
point  2  *  +1.  This  situation  is  shown  in  Figure  1,  where 
the  three  roots  are  designated  z\,  Z2,  and  23  =  "Z2 
(the  overbar  is  used  to  indicate  the  complex  conjugate  of  a 
quantity) : 


zx  ~  d, 

and 

z2  =  pe1^  *  x  +  i  y, 

where  d,  x,  and  y  are  real  quantities  and 

,  2  .  2.  J5 

p  «  (x  +  y  )  , 

<)>  =  tan-1  (y/x)  . 


(1) 

l 

I 


(2) 


(3) 

(4) 
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r'fe,.- 


The  complex  variable,  z,  is  defined  in  the  usual  way. 


z  —  r  e  ,  (5) 

where  r  is  the  modulus  and  the  argument  6  is 
defined  as 


e 


COS  1  (a)  T 
n 


(6) 


The  symbol  co  represents  the  natural  frequency  associated 
with  a  complex  pair  of  roots,  T  represents  the  sampling 
period,  and  c  represents  the  damping  ratio  associated  with 
the  pair  of  complex  roots.  When  0  equals  0  or  it,  then  z  is 
a  real  (rather  than  complex)  variable. 

If  it  is  assumed  that  z\,  Z2,  and  Z3  comprise  the 
roots  of  a  particular  system,  then  the  characteristic 
equation  (C.E.)  may  be  written  as 


C.E.  (z)  =  (z-z^  (z-z2)  (z-z3) 

=  (z-d)  (z-pei(^)  (z-pe  i<*>) 

=  z3  -  (d  +  2  p  cos  <J>)  z2  +  (p2  +  2d  p  cos  4>)  z 
2 


-  dp‘ 

=  z3  -  (d  +  2x)  z2  +  (p2  +  2dx)  z  -  dp2  =  0 


(7) 


One  observes  that  as  the  locations  of  the  three  roots 
approach  each  other,  i.e.,  as  p  approaches  d,  one  obtains 
(in  the  limit)  from  Equation  (7)  the  expression 


lim  C.E.  (z)  =  lim  C.E.  (z)  =  (z-d)3  =  0,  (8) 

p-*-d  x-*-d 

y-+0 


as  expected.  Computer  and  scaling  problems  can  occur  if  d 
approaches  the  absolute  value  of  unity,  as  is  explained 
below. 
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As  can  be  seen  with  the  aid  of  Figure  1,  when  d 
approaches  either  +1  or  -1,  the  three  (in  this  example) 
clustered  roots  approach  one  of  the  two  singularities 
(i.e.,  z  *  +1,  z  =  -1,)  in  the  z-domain.  If  it  is  desired 
to  solve  the  system  characteristic  equation  for  its  roots 
(the  so-called  eigenvalue  problem)  in  order  to  assess 
system  stability,  computational  techniques  often  run  into 
difficulty  if  the  roots  (eigenvalues)  are  close  to  each 
other.  Sometimes  this  problem  can  be  circumvented  by 
increasing  the  scale  of  the  system,  i.e.,  -  in  the  case  of 
Figure  1  -  by  increasing  the  scale  of  the  figure.  However, 
because  the  origin  of  the  vectors  to  the  roots  is  at  z  =  0 , 
when  the  roots  are  clustered  closely  and  on  the  unit  circle 
(such  as  at  z  *  +1),  it  is  difficult  to  separate  them  suf¬ 
ficiently  for  the  computer  to  solve  the  associated 
eigenvalue  problem.  It  can  be  induced  that  if  the  problem 
at  hand  has  even  more  (than  three)  clustered  roots,  the 
computational  problem  becomes  more  severe.  Thus  a  means  of 
avoiding  this  "pathological"  (computer  jargon)  problem  is 
to  be  desired  if  it  appears  that  such  problems  might  arise 
in  actual  practice. 

IMAGINARY 


REAL 


Figure  1.  Clustered  roots  in  z-domain. 
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B.  Likelihood  of  Occurrence 


Before  attempting  to  attack  the  problem  of  clustered 
roots,  it  is  desirable  to  assess  the  possibility  of  their 
actually  occurring  in  practice.  Digital  design  experience 
shows  that  the  probability  of  system  roots  occurring  at  the 
z  =  -1  singularity  is  unlikely,  unless  deliberately  placed 
there  by  a  digital  filter.  Since  the  system  designer  has 
the  option  of  removing  such  roots  by  altering  the 
characteristics  of  his  digital  filter,  this  problem  can  be 
solved  readily.  However,  the  problem  is  not  so  readily 
solved  if  the  roots  are  clustered  at  the  z  =  +1  singular¬ 
ity. 


Suppose  that  a  system  plant  has  a  second  order 
characteristic,  such  as  one  would  find  typically  with  a 
bending  mode  or  in  the  denominator  of  a  transfer  function 
(T.F.)  describing  the  characteristics  of  some  sensors  or 
actuators.  A  typical  transfer  function  might  be 


T.F.  (s)  =  N(s) _  ,  (9) 

(s  +  a)2  +  w2 

where  the  numerator  N(s)  may  have  any  form  (in  this  discus¬ 
sion).  An  example  for  a  bending  mode  would  have  N(s) 
represent  the  modal  gain,  a  would  be  the  modal  damping 
ratio,  and  u  would  be  the  square  root  of  the  difference  of 
the  squares  of  the  modal  natural  frequency  (wn)  and  a, 
i.e. 


U>  =  U2  -  a2  )%  .  (10) 

The  z-transform  of  the  transfer  function  of  Equation  (9)  is 
seen  to  be 


T.F.  (z) 


_ Ms) _ 

z2  -  2ze~aTcos  u)T  +  e"2aT 


(ID 


where  T  represents  the  sampling  period  of  the  digital 
system.  The  nature  of  numerator  N(z)  depends  on  the  form 
of  N(s).  The  important  thing  to  observe  is  that  if  a  or  T 
or  both,  (i.e.,  if  their  product)  approach  zero,  then  the 


denominator  of  Equation  (11)  can  approach  the  form. 


D(z)  =  (z-1)2. 


(12) 


In  this  unfortunate  case,  the  roots  would  be  clustered  near 
(or  on,  in  the  limit)  the  singularity  at  z  =  +1,  giving 
rise  to  the  problems  brought  forth  in  Section  A  above. 

Thus  it  unfortunately  appears  possible  for  the  clustering 
problem  to  actually  arise  in  a  physical  system.  The 
probability  of  such  an  occurrence  depends  on  the  actual 
dynamic  nature  of  the  plant's  sensors,  actuators,  and 
dynamic  bending  properties.  Judicious  engineering  design 
can  of  course  obviate  such  problems  brought  forward  because 
of  actuators  or  sensors,  if  desired.  However,  changing  the 
bending  characteristics  of  the  missile  plant  can  be 
difficult.  All  is  not  lost,  however,  because  it  is  the 
product  of  a  and  T  in  Equation  (11)  that  gives  rise  to  the 
problem.  Hence,  again  the  designer  can  select  the 
numerical  value  of  T  that  will  preclude  such  a  clustering 
problem.  If  other  design  constraints  conspire  against  the 
digital  system  designer,  he  might  desire  another  analytical 
approach  that  might  solve  the  clustering  problem.  Should 
that  desire  materialize,  another  approach  is  proposed  and 
discussed  in  the  following  sections  of  this  report. 

III.  PROPOSED  SOLUTION 

The  problem  investigated  in  Section  II  is  attacked  in 
this  section  by  applying  a  transformation  from  the  complex 
z-domain  to  the  complex  w-domain.  As  desired  in  Paragraph 
2.2  of  the  referenced  Technical  Requirement,  it  is  intended 
that  this  technique  "mathematically  spread  the  system  poles 
and  zero  via  intermediate  complex  transform  techniques  to 
eliminate  the  sensitivity". 

As  in  Section  II  (for  the  sake  of  simplicity),  one  may 
again  assume  the  dynamics  of  a  system  are  represented  by 
the  three  roots  of  Equations  (1)  -  (2).  The  usual 
definition  used  in  mapping  from  the  z-domain  to  the 
w-domain  is: 


Z  =  11  i  w). 

(1  -  w) 


(13) 
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The  complex  variable  w  is  composed  of  real  and  imaginary 
parts  which  may  be  denoted  in  the  usual  manner. 


w  =  a  +  ito  . 
w  w 


Applying  Equation  (13)  to  root  t.\  (1),  one  may  obtain  its 
equivalent  location  in  the  w-domain: 


_  (1  +  w)  _  _ 

-  VI - —ITT  —  Z,  —  d. 


If  one  solves  Equation  (13)  for  w,  one  obtains  the  root 
locat ion 


...  _  _(1  -  d) 

W1  “  71+  d) 


It  may  be  observed  in  passing  that  the  real  root  location 
in  the  z-domain  has  been  transformed  into  a  real  root  loca¬ 
tion  in  the  w-domain.  In  a  similar  manner  the  equivalent 
w-domain  expression  for  the  complex  conjugate  root  of  Equa¬ 
tion  (2)  may  be  obtained: 


2  =  (1  -  w)  -  *2  =  C  e1*  =  *  +  iy 


When  rearranged  and  solved  for  w,  one  may  obtain 


w2  =  - 


(1  -pei<p) 
(1  +pe1(^) 


=  -U  -  x  -  y  )  + 
(1  +  x)'  +  y 


)  +  i  2) 
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The  expression  for  W3  is  the  complex  conjugate  of  Equa¬ 
tion  (18).  The  three  roots  may  now  be  mapped  onto  the  com¬ 
plex  w-domain  as  shown  in  Figure  2. 

If  it  is  assumed  that  w^,  W2f  and  W3  comprise  the 
roots  of  the  same  system  analyzed  in  Section  II,  then  the 
characteristic  equation  may  be  written  as 


C.E.  (w)  =  (w-w^  (w-w2)  (w-w3)  =  0.  (19) 

Alternatively,  the  characteristic  equation  may  be 
determined  by  subsituting  the  value  for  z  of  Equation  (13) 


IMAGINARY 


REAL 


Figure  2.  Clustered  roots  in  w-domain. 
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into  C.E.(z)  of  Equation  (7)  to  obtain 


C.E. (w)  = 
+ 


+ 


<l+d+2x+  p2+2dx+d  p2)w3 
(3+d+2x-  p2-2dx-3d  p2)w2  + 

(l-d-2x+  p2+2dx-d  p2)  =  0. 


(3-d-2x-  p2-2dx+3d  p2)w 

(20) 


One  observes  in  this  case  that,  as  the  locations  of  the 
three  roots  coalesce,  i.e.,  as  p  approaches  d,  one  obtains 
(in  the  limit)  from  Equation  (20)  the  expression 

lim  C.E.  (w)  =  lim  C.E. (w)  =  (1  +  d) 3w3  +  3(1  +  d)2(l  -  d)w2 
p+d 

x-»-d  +  3  (1+d)  (1-d)  2w  +  (1-d)  3  =  0. 

Y+°  (21) 


If  d  approaches  unity  (as  before),  then 
lim  C.E. (w)  =  (2w) 3  =  0, 

d-1  (22) 


and  the  three  roots  approach  (as  a  cluster)  the  origin  of 
the  w-plane.  Now,  however,  if  the  scale  of  the  system  - 
in  the  case  of  Figure  2  -  is  increased,  success  is  met  in 
spreading  apart  the  roots  because  the  vector  to  each 
originates  at  the  nearby  origin.  Recall,  in  the  z-domain, 
increasing  the  scale  did  not  substantially  spread  the  roots 
apart  because  the  vectors  to  the  roots  went  (approximately) 
from  the  origin  to  the  z  =  +1  point. 

IV.  COMPARISON 

The  scaling  problem  and  proposed  solution  alluded  to 
above  in  the  w-domain  may  not  be  immediately  apparent. 
Therefore  several  numerical  comparisons  are  made  in  this 
section.  Ten  runs  were  made  for  selected  values  of  roots 
in  the  z-domain.  The  corresponding  roots  in  the  w- 
domain  were  calculated  for  comparison.  The  results  were 
obtained  using  the  program  of  the  Appendix  on  the  Hewlett- 
Packard  9100B  Calculator.  They  are  summarized  in  Table  1. 


TABLE  1. 


ROOT  LOCATIONS  IN  z-  and  W-DOMAINS 


Run 

Z!*^ 

z2*x+iy 

Wi 

w2 

1 

0.9 

0 . 8+ iO . 2 

-0.05 

-0 . 10+i0 . 12 

2 

0.96 

0.95+i0 .02 

-0.02 

-0.01+i0.03 

3 

0.98 

0.97+i0,02 

-0.010 

-0.015+i0.010 

4 

0.99 

0 .98+iO.Ol 

-0.005 

-0.010+i0.005 

5 

0.98 

0.99+i0.01 

-0.010 

-0.005+i0.005 

6 

0.996 

0.95+0.001 

-0.002 

-0.00251+ i5X10“4 

7 

0.995 

0 .996+iO .001 

-0.00251 

-0 . 00200+ i5X10“4 

8 

0.999 

0 .999+iO .001 

-5X10-4 

-5X10-4+ i5XlO“4 

9 

0.9999 

0.9999+i0.0001 

-5X10-5 

-5X10“5+i5X10-5 

10 

0.99999 

0.99999+0.00001 

-5X10-5 

-5X10-5+15X10-5 

It  is  seen  that  the  distance  between  the  z  roots  and  the 
singular  point  at  z  ■  +1  is  about  twice  the  distance 
between  the  w  roots  and  the  origin  of  the  w-plane,  even 
though  (in  each  case)  the  distances  are  quite  small.  How- 
ever,  the  effect  of  increasing  the  scaling  in  the  z-domain 
is  not  nearly  so  great  as  the  effect  of  increasing  the 
scaling  in  the  w-domain.  This  is  seen  by  example  on  Figure 
3  where  the  z-  and  w-  domain  roots  are  plotted  on  an 
Tncreases  scale. 


IMAGINARY 


r-DOMAIN 


REAL 


X 


W- DOMAIN 


IMAGINARY 
1 iO.OOS 


— i - 

-0.010 


X 


-0.005 


REAL 


X 


10.005 


Figure  3.  Clustered  roots  of  run  #4. 

(Shows  effect  of  scaling) 


It  is  readily  observed  that  increases  in  scaling  in  the 
z-domain  do  not  have  a  significant  effect  on  roots 
clustered  about  the  z  =  +1  singularity.  This  is  because 
the  vectors  to  these  roots  emanate  from  origin.  The  same 
is  not  true  in  the  w-domain.  Here  it  is  seen  that  large 
increases  in  scale  have  profound  effect  on  spreading  out 
the  roots  clustered  about  the  origin.  This  is  because  (in 
this  case)  the  vectors  to  these  roots  emanate  from  the 
origin  to  the  (in  this  case)  nearby  root  locations. 
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V.  CONCLUSION 

An  investigation  has  been  marie  to  investigate  the 
potential  existence  of  a  clustering  of  system  poles  and 
zeros  in  the  z-domain.  It  was  shown  by  example  in  Section 
IIA  how  such  a  problem  might  arise.  In  Sectib^IIB  the 
likelihood  of  such  an  occurrence  was  assessed .— J It  was 
shown  that  such  a  situation  could  occur,  but  that  there  are 
design  techniques  to  prevent  it.  If  these  fail,  an 
alternate  approach  is  described  4n- Sect  ion  whereby  the 
system  poles  and  zeros  are  transformed  from  the  complex 
z-domain  to  the  complex  w-domain.  It  then  is  possible  to 
increase  the  scaling  (in  the  w-domain)  with  significant 
effect  in  separating  the  clustered  poles  and  zeros,  thereby 
obviating  computational  difficulties  that  might  have  arisen 
because  of  the  clustering  in  the  z-domain. 
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